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Activation of skeletal muscle fibers requires rapid sarcolemmal action potential (AP) conduction to ensure uni- 
form excitation along the fiber length, as well as successful tubular excitation to initiate excitation-contraction 
coupling. In our companion paper in this issue, Pedersen et al. (2011./. Gen. Physiol. doi:10.1085/jgp. 201010510) 
quantify, for subthreshold stimuli, the influence upon both surface conduction velocity and tubular (t) -system ex- 
citation of the large changes in resting membrane conductance ( Gm) that occur during repetitive AP firing. The 
present work extends the analysis by developing a multi-compartment modification of the charge-difference 
model of Fraser and Huang to provide a quantitative description of the conduction velocity of actively propagated 
APs; the influence of voltage-gated ion channels within the t-system; the influence of t-system APs on ionic homeo- 
stasis within the t-system; the influence of t-system ion concentration changes on membrane potentials; and the 
influence of Phase I and Phase II G M changes on these relationships. Passive conduction properties of the novel 
model agreed with established linear circuit analysis and previous experimental results, while key simulations of AP 
firing were tested against focused experimental microelectrode measurements of membrane potential. This study 
thereby first quantified the effects of the t-system luminal resistance and voltage-gated Na + channel density on sur- 
face AP propagation and the resultant electrical response of the t-system. Second, it demonstrated the influence of 
G M changes during repetitive AP firing upon surface and t-system excitability. Third, it showed that significant K + 
accumulation occurs within the t-system during repetitive AP firing and produces a baseline depolarization of the 
surface membrane potential. Finally, it indicated that G M changes during repetitive AP firing significantly influence 
both t-system K + accumulation and its influence on the resting membrane potential. Thus, the present study 
emerges with a quantitative description of the changes in membrane potential, excitability, and t-system ionic ho- 
meostasis that occur during repetitive AP firing in skeletal muscle. 



INTRODUCTION 

Recent studies have demonstrated that the resting mem- 
brane conductance (G M ) of skeletal muscle is highly 
regulated during repetitive firing of short trains of ac- 
tion potentials (APs) that replicate excitation patterns 
occurring during activity in skeletal muscle (Pedersen 
et al., 2009a,b) . In fast-twitch rat muscle, this regulation 
of G M has two distinct phases. At the onset of AP firing, 
Phase I involves the inhibition of C1C-1 channels through 
a PKC-mediated mechanism causing a decline in G M to 
^40% of its resting value. Then, during prolonged ac- 
tivity, Phase II changes involve the opening of C1C-1 and 
K ATP channels, increasing G M to four to five times its 
value in quiescent fibers. Upon cessation of AP firing, 
G M recovers to its level before AP firing in 1-5 min. 

In our companion paper (see Pedersen et al. in this 
issue) , we demonstrated the significance of this G M reg- 
ulation for subthreshold electrical properties in muscle 
fibers, and from this predicted the effects of such G M 
changes upon their excitability. In particular, our study 
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described the relative sensitivity of the different aspects 
of muscle excitability, including neuromuscular trans- 
mission, sarcolemmal AP propagation, and tubular 
(t)-system excitation, to such regulation. The study used 
a linear circuit analysis appropriate for subthreshold 
electrical membrane phenomena. This approach first 
allowed for the development of analytical solutions for 
three cable models of muscle fibers. To determine which 
of these models gave the best representation of the elec- 
trical properties of rat extensor digitorum longus (EDL) 
muscle fibers in which the G M regulation has been ob- 
served, experimental measurements of membrane im- 
pedance properties were compared with the electrical 
characteristics of the cable models. This demonstrated 
that circuit models of rat EDL muscle fibers require a 
substantial luminal t-system resistance to account for 
experimental observations of their impedance prop- 
erties and the velocity with which APs propagate in 
these fibers. 
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The study went on to predict that such a luminal resis- 
tance enhances the propagation velocity of sarcolemmal 
APs. It also showed that a luminal resistance in series 
with the t-system membrane divides the voltage gradi- 
ent between the intracellular and interstitial spaces into 
voltage drops across the t-system membrane and the lu- 
minal resistance. Because this voltage division is highly 
frequency dependent, the luminal resistance has impor- 
tant consequences for t-system excitation. Thus, high- 
frequency current will predominantly generate voltage 
gradients across the luminal resistance, whereas lower 
frequency components will be important for t-system 
excitation. The analytic approach of the previous study 
therefore demonstrated that Gm changes predominandy 
affect the low-frequency membrane impedance. It was 
further demonstrated, by convolving experimental APs 
with the circuit models, that t-system excitation is a low- 
frequency phenomenon that is highly dependent on 
Gm- Collectively, our companion study thus demon- 
strates that models of rat EDL muscle fibers must include 
a luminal resistance and suggests that neuromuscular 
transmission and t-system excitation are more sensitive 
to Gm regulation than is sarcolemmal AP propagation. 

The linear circuit analysis used in our companion 
study was also useful in providing analytic expressions 
that could distinguish the appropriate equivalent cir- 
cuit model for rat EDL muscle fibers and quantify the 
t-system luminal resistance. However, such a linear analy- 
sis is necessarily incomplete when exploring the possi- 
ble physiological roles of the t-system luminal resistance 
and G M regulation for the nonlinear membrane phe- 
nomena involved in AP propagation. Such nonlinear 
properties, explored in the present study, include volt- 
age-dependent, time-dependent, and rectifying and 
Na + /K + -ATPase-mediated currents, as well as alterations 
in intracellular and t-system ionic concentrations dur- 
ing repetitive activity. 

The aim of the present study was therefore to quan- 
tify the influence of t-system luminal resistances and Gm 
regulation in skeletal muscle excitability and t-system 
ionic homeostasis. It uses a nonlinear, iterative approach 
based upon the charge-difference (CD) model of Fraser 
and Huang (2004) and Fraser et al. (2005a,b). The de- 
velopment of a new model was necessary because no ex- 
isting model simulates the full range of physical and 
electrophysiological properties that underlie t-system 
ionic homeostasis and its relationship with the mem- 
brane potential. Thus, early models applying circuit 
theory with voltage- and time-dependent conductances 
allow for the simulation of individual APs in the absence 
of ionic concentration or osmotic changes and are there- 
fore unsuitable for simulating trains of APs where ionic 
concentrations shift substantially (Adrian and Peachey, 
1973; Adrian and Marshall, 1976). Subsequent, more 
realistic models permitting the simulation of APs within 
the restricted extracellular space of a whole muscle 



(Henneberg and Roberge, 1997), and of t-system K + 
handling during AP firing (Wallinga et al., 1999), do 
not incorporate osmotic water movements. Therefore, 
they do not reach unique steady-state solutions that are 
independent of the initial values of key modeled vari- 
ables, such as intracellular and intra-t-system ion con- 
centrations (Fraser and Huang, 2007) . 

The modeling approach adopted here encompasses 
the known determinants of AP propagation and reaches 
a true history-independent steady state, thereby remain- 
ing applicable despite large perturbations in ionic con- 
centration. It incorporates terms describing voltage- and 
time-dependent ion conductances and Na + /K + pump 
activity; represents muscle fiber surface geometry using 
99 linearly connected fiber segments permitting simula- 
tion of surface conduction; represents tubular geome- 
try in terms of 20 concentric shells per fiber segment, 
each separated by a small luminal series resistance as 
was experimentally verified in our companion study 
(Pedersen et al., 2011); and simulates osmotic water 
fluxes, thereby permitting it to reach a unique history- 
independent steady state (Fraser and Huang, 2007). 

History independence of the modeled variables was 
demonstrated by initiating the new model from several 
sets of widely divergent unphysiological values for all 
ion concentration variables. In each case, the model re- 
laxed to an identical and physiologically reasonable 
steady state, thereby demonstrating that it was capable 
of investigating the determinants of intracellular and 
intra-t-system ionic homeostasis, and additionally allow- 
ing for the use of model-derived initial values for vari- 
ables such as intra-t-system ion concentrations for which 
there is little available experimental data. 

The initial application of this model confirmed the 
prediction from the preceding analytical study that the 
t-system luminal resistance enhances the sarcolemmal 
AP propagation velocity and reduces t-system excita- 
tion. It further demonstrated the influence of active AP 
generation within the t-system on t-system excitation 
and t-system ionic homeostasis, and the effect of these 
upon the surface membrane potential. It was then used 
to explore the role of Gm regulation in determining the 
threshold stimulus for AP firing, for t-system excitability, 
and for the homeostasis of K + ions within the t-system 
during repetitive AP firing. The outputs of simulated 
trains of APs were shown to be in close quantitative 
agreement with experimental intracellular recordings 
of AP trains obtained from rat EDL muscle fibers. 

Theory 

Development of a CD model of rat skeletal muscle 

The electrical properties of rat EDL skeletal muscle fi- 
bers were modeled on CD principles (Fraser and Huang, 
2007), adapting and extending the model of Fraser and 
Huang (2004) and Fraser et al. (2005a,b) to include key 



96 Membrane conductivity and muscle excitability 



nonlinear membrane phenomena of AP excitation and 
conduction and t-system excitation, as depicted in Fig. 1. 
The new model thus simulated the following features: 
longitudinal subdivision of the model fibers into multi- 
ple subsections connected to form a cable model; subdi- 
vision of the t-system into a series of concentric shells 
separated by series resistances representing the t-system 
luminal resistance; voltage- and time-dependent ion chan- 
nels within the surface and t-tubular membranes; and 
ion and osmotic water fluxes across the surface and 
t-system membranes and the t-system luminal resistances. 

The model parameters were calibrated to the passive 
electrical properties of rat fast-twitch EDL muscle fibers 
that were obtained using the analytic model in our com- 
panion study (Pedersen et al., 2011). This allowed its 
subsequent use in evaluating the physiological signifi- 
cance of the t-system luminal resistance and Gm regula- 
tion in active muscle fibers that has been observed 
previously in AP-firing EDL fibers (Pedersen et al., 2009a) . 
The model was developed to adhere stricdy to four key 



principles: (a) current continuity; (b) conservation 
of charge and mass; (c) electrodiffusion of ions; and 
(d) osmotic equilibrium at steady state. 

(a) The principle of current continuity across the resistances 
and capacitances of the CD model. The simulation of mul- 
tiple apposing finite compartments within strict conser- 
vation principles requires that any flow of current 
completes a circuit to ground, and also that such cur- 
rents, whether capacitative or resistive, correlate pre- 
cisely with the correct changes in compartmental ion 
contents. Both current continuity and conservation 
principles can be achieved straightforwardly for ionic 
fluxes through ion channels and across resistances by 
the following equations: 

A[«L = /sAn, 0 )A*, (moirV 1 ), 
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Figure 1. The ionic fluxes and capacitances simulated in the CD model of skeletal muscle cable properties. A multi-compartment 
model was developed to model the cable properties of skeletal muscle. The muscle fiber was divided into 99 longitudinal segments, the 
length of which could be varied. Each cable segment contained a t-system compartment that could be simulated as a single compartment 
or further subdivided into several concentric shells. In each cable segment, ionic fluxes through background and voltage-gated ion chan- 
nels, calculated using the Goldman equation, and Na + /K + -ATPase (pump) fluxes were simulated across the sarcolemma membrane (/ s ), 
and across the membrane of each t-system shell (/„). Ionic fluxes across the t-system access resistance (Jt^o)), between t-system shells 
(/( n _>n+i))> an d between adjacent cable segments of the muscle fiber (/ (1) ) were calculated using an electrodiffusion equation according to 
the prevailing concentration and electrical gradients. In addition to the depicted ionic currents, water fluxes were also modeled, allow- 
ing the model to reach a steady state that is independent of the initial concentration of any ion. Potential differences were calculated, as 
described in the Theory section, across the surface membrane capacitance (C m ), across the membrane capacitance of each t-system shell 
(C„), between each cable segment, between the extracellular space and the outer shell of the t-system, and between each t-system shell. 
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where [S] represents the concentration of any solute 
(mol 1 ), the subscripts a and b denote the compart- 
ments between which the ionic flux, J (mol cm 2 ), is 
flowing, and A m is the membrane area per unit volume 
(cm 2 1~ ) for each compartment. Conservation of charge 
and ion concentration is clear from the relationship 
A[S]yA m(a) = -A[S] s /A m(s) . 

The CD modeling approach additionally permitted 
calculation of the less straightforward requirement to 
characterize the precise ionic concentration changes re- 
sulting from capacitive currents, in which association of 
one ion species with one "plate" of a capacitor might 
displace or attract different ion species from or to the 
other plate, such that A.[S\ a /A m{a) * — A[S] b / A mW , in 
marked contrast to the situation with currents through 
ion channels. For example, hyperpolarization of the 
t-tubular membrane could result in the association of 
Na + and/ or CI dissociation with the outer membrane 
surface, but the dissociation of K + from the inner mem- 
brane surface. This possibility is made explicit by the 
following relationships for the surface and t-tubular 
membrane capacitance (C m , F cm 2 ) in a muscle fiber 
with a homogenous t-tubular compartment separated 
from the bulk extracellular fluid (ECF) by an access re- 
sistance, Ra: 

E m =Q sm /{C m A m ), (V) 

and: 

£ t = Q m /(c m A)> (v) 

while: 

E A = E m -E t , (V) 

where Q m is the net charge associated with the surface 
membrane capacitance, Q m is the net charge associated 
with the t-system membrane capacitance, with each ex- 
pressed per liter of cell volume (C P 1 ), and A m and A t 
represent the membrane areas of the sarcolemma 
and t-system, respectively, referred to total cell volume 
(cm 2 . £m represents the trans-sarcolemmal membrane 
potential, Et is the trans-tubular membrane potential, 
and Ea is the potential across the tubular access resistance. 
Ionic fluxes then direcdy influence the ion concentrations 
and hence the net charges within the intracellular (i) and 
t-tubular (t) compartments. These relate to the charges 
associated with membrane capacitances by: 

Qi = Q m + a» ! (cr 1 ) 
&= -^-(ci" 1 ) 

These relationships encapsulate the principles that 
(a) intracellular charge may be associated with either 



Cm or Q; and (b) any deposition of charge on one plate 
of Q must be matched by an equal charge flowing from 
its opposing plate. These relationships incorporate the 
assumptions that (a) the outer surface of C m is grounded 
via a pathway of insignificant resistance; and (b) the outer 
surface of Q is grounded via a nonzero access resistance, 
Ra, thereby allowing and Et to differ. Collectively, 
these assumptions imply that any alterations in Q, for ex- 
ample because of currents across Ra, are matched by equal 
alterations in charge on Qn, thereby completing a cir- 
cuit to ground. 

(b) Strict conservation principles in CD modeling of multiple 
interconnected finite compartments. By tracking concen- 
tration changes and then calculating electrical poten- 
tials directly from the concentrations of charge carriers 
and the associated capacitance terms, the CD approach 
results in a model that obeys fundamental charge and 
mass conservation principles. This allows the model to 
reach true steady states that are independent of the ini- 
tial values of any modeled variable (Fraser and Huang, 
2007). Thus, for a single membrane-bound compart- 
ment, the electrical potential with respect to the bulk 
ECF is given at any time by: 

E m = Qi/C m , (V) 

where is the membrane potential, and C m is the 
membrane capacitance. Q_ denotes net charge, and the 
subscript "i" denotes the intracellular space, such that: 

Q i= F([Na + ] 1 +[K + ] 1 -[Cr] 1 + z x [4), (CI" 1 ) 

where z x is the mean charge valency of the various mem- 
brane-impermeant solutes, denoted X, and Fis Faraday's 
constant. The influence of z x and [X] ; on the steady state 
of the model has been described previously (Fraser and 
Huang 2004, 2007; Fraser et al., 2005a) . 

With the current-continuity assumptions explored 
above, relationships may be given for and E t in terms 
of the CDs within each compartment. First, if the t-system 
is considered as a single homogenous compartment, par- 
tially bounded by a membrane with capacitance per 
unit area that is identical to that of the surface mem- 
brane, then: 

F _fw) m 

K ~ (cjl) ' (V) 

E c =^, (V) 

[Ql = F([Na\ + [K\ - [Cl\ + z x [X\) 

[Q\ = F([Na\ + [K\ - [CJ\ + [Z\ ) , CI" 1 
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where V t is the volume of the t-system relative to the cell 
volume, V c ; and Z (t) denotes the fixed charge density 
within the t-tubules (C 1~ ). It is implicit in these equa- 
tions that the ionic activities of the mobile ions are 1. 
Sign conventions are as follows: the surface membrane 
potential (E m ) is expressed as the potential of the intra- 
cellular space with respect to the ECF; and the t-system 
membrane potential (£i( n )) is expressed as the potential 
of the intracellular space with respect to the potential of 
the t-system lumen in each shell, n. 

To calculate transmembrane potentials when the 
t-system is modeled as multiple concentric shells, the ex- 
cess charge within the t-system is simply summed, thus: 



-[Q] tM y» 



(V) 



where V n (dimensionless) is the volume of each t-system 
shell, n, relative to the cell volume, and V t is the volume 
of the t-system relative to cell volume. 

(c) Electrodiffusion of ions. In all cases, ionic movements 
were calculated according to both the concentration 
gradient and the local electrical gradient. This contrasts 
with previous models of skeletal muscle fibers that con- 
sider either diffusional fluxes or electrical currents across 
t-system luminal resistances, but not both (Falk and 
Fatt, 1964; Freygang et al., 1964; Adrian and Peachey, 
1973; Adrian and Marshall, 1976; Henneberg and 
Roberge, 1997). However, preliminary simulations dem- 
onstrated that AP activity resulted in alterations in both 
t-system ion concentrations and t-system electrical po- 
tential, thereby necessitating the simulation of both dif- 
fusional and electrically driven ion movements. 

Transmembrane ionic fluxes were calculated using 
the Goldman flux equation (Goldman, 1943). Fluxes 
across R A and between t-system shells were calculated as 
the sum of electrical drift and diffusional fluxes: 



• + f([S] ( „-i) -[%,). (mol s" 1 ) 



'[C]„ 



where cr t is the t-system tortuosity factor (Wallinga et al., 
1999) (dimensionless) ; A n is the area of shell n (dm 2 ) ; A 
is the distance between successive t-system shells (dm) ; 
g t is the inter-shell conductance (S dm 1 ) for shell n; 
[S] n /[C] „ is the concentration of S relative to the total 
mobile ion concentration within shell n (dimensionless) , 
such that each ionic species, S, contributes an appropriate 



proportion of the total conductance; and D is the diffu- 
sion coefficient (dm 2 s _1 ), assumed to be equal for all 
modeled ions. 
Fluxes across Ra were calculated as follows: 



/S(e^O) 



■An 
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mol! 



where the D/ Rp^g t term ensures that the ratio between 
the diffusive and electrical components of the flux 
across Ra is equal to that between t-system shells. 

(d) Osmotic equilibrium at steady state. Fraser and Huang 
(2004, 2007) showed that when the steady state of a 
CD model cell is constrained by the requirement 
for osmotic equilibrium at steady state, then the steady- 
state values of these variables were uniquely defined 
by the permeabilities and extracellular concentrations 
of each ion and by the Na + /K + pump density. Further- 
more, certain concentrations are tightly constrained 
at osmotic equilibrium. For example, when CP is the 
only extracellular anion and therefore contributes 
half of the extracellular osmolality, [Cl~] t must be 
equal to [CP] e for t-system osmolality to be equal to 
extracellular osmolality. Thus, osmotic water move- 
ments are fundamental to ensure accurate steady- 
state behavior that is constrained to what is physically 
possible. This allows CD modeling to be used to simu- 
late manipulations expected to have a large influence 
on ion concentrations, in contrast to models that 
do not reach true steady states (Fraser and Huang, 
2007). Thus, transmembrane water fluxes were calcu- 
lated as: 



h 



H20(a->b) 



where 17 is the osmolality of each compartment and 
Pmo is the hydraulic permeability of the membrane, 
taken as 156 pm s (Frigeri et al., 2004). Water fluxes 
across R A and between t-system shells were assumed to 
be instantaneous to maintain osmotic equality in all ex- 
tracellular compartments. 

It is clear from the profound changes in t-system vol- 
ume with manipulations of extracellular ion concen- 
trations or osmolality (Rapoportetal., 1969; Dulhunty, 
1982; Chawla et al., 2002) that the t-system cannot be 
treated as a rigid structure. The model therefore 
allowed the t-system volume to be determined by the 
balance of osmotic fluxes across the t-system mem- 
brane and Ra- This gave stable steady-state values when 
a small fixed negative charge term was included in the 
t-system and responded appropriately to changes in 
extracellular osmolality in preliminary experiments 
(unpublished data) . 
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Simulation of cable properties and AP firing 

A length of muscle fiber was simulated as 99 identical 
segments. Ion concentrations and membrane potentials 
were calculated for each segment as described above. 
Inter-compartment ion fluxes were calculated for each 
intracellular ionic species according to the prevailing 
electrical gradient between segments: 

/ S ™,=^r[ Tl J.^oidm s ) 

where A x is the membrane area of compartment x (dm 2 ) , 
V x is its volume (dm 3 ), R L is the longitudinal resistance, 
and AL is the distance between the centers of adjacent 
compartments. As each segment was identical, ion con- 
centrations did not differ significandy between adjacent 
segments; therefore, it was not necessary to include a 
diffusional flux term. 

The simulations were performed using a 100-pm seg- 
ment length. Initial simulations to determine the optimal 
segment length compared the behavior of model fibers 
with 100-pm segment lengths to those with 10-pm seg- 
ment lengths. Over a range of frequencies from 25 to 
3,000 Hz, frequency-velocity relationships for passive si- 
nusoidal currents were not influenced by this change in 
segment length (difference in velocity <\% at each fre- 
quency; not depicted) . 



APs were modeled by adding voltage-gated Na + 
and K + channels to the surface and t-tubular mem- 
branes. Voltage-gated ion channels were simulated 
using Hodgkin-Huxley gating equations (Hodgkin and 
Huxley, 1952), with conductances and gating parame- 
ters calibrated to rat skeletal muscle (Cannon et al., 
1993; Henneberg and Roberge, 1997; Wallinga et al., 
1999). Nevertheless, the present model differs signifi- 
cantly from these existing models in several important 
respects. First, all modeled variables, including t-system 
ion concentrations and ion channel-gating variables, 
reach stable steady-state values that are independent 
of initial conditions. In contrast, these existing models 
treat some or all ion concentrations as fixed parame- 
ters. Second, the present model simulates ion diffu- 
sion between the t-system and extracellular space 
according to the prevailing electrical and concentra- 
tion differences. In contrast, existing models simulate 
passive diffusion, yet particularly when E t lags E m 
during the early phase of the surface AP, the voltage 
between the t-system and the ECF was shown to reach 
almost 100 mV in these earlier studies (Wallinga 
etal., 1999). 

A stimulus was required to initiate APs. This was ap- 
plied to the central compartment in the cable model. 
To maintain strict conservation principles, stimulation 
currents were simulated as fluxes of K + and CP of equal 



TABLE I 

Summary of physical and electrophysiological parameters used 
Symbol Description Value Reference 

Physical parameters 

A m Surface membrane area 769231 cm 2 1~ 1 



A, 


t-System membrane area 


3076924 cm 2 l" 1 


1 


i; 


Relative t-system volume 


3.2e-3 liter.r 1 


2,7 




Membrane capacitance 


le-6 F cm -2 


1 


K 


t-System access resistance 


72.3 a cm 2 


1 


n, 


Longitudinal fiber resistance 


180 0 cm 


1 




t-System luminal conductance 


3.7e-3 S cm" 1 


2 


ft 


t-System tortuosity factor 


0.21 


2 


Electrophysiological parameters 








N, 


Na + /K + -ATPase sarcolemmal density 


9e-12 cnr 2 


3 


At 


Na + /K + -ATPase tubular membrane density 


4.5e-12 cm~ 2 


3 


fka 


Background Na + permeability 


2.16e-8 cm s" 1 


4, 5 


ft 


Background K + permeability 


6.3e-7 cm s _1 


4, 6 


Pa 


Background CI" permeability 


2.7e-6 cm s~ 


4, 6 




Voltage-gated Na + channel sarcolemmal permeability 


0.00265 cm s" 1 


7 




Voltage-gated Na* channel t-system permeability 


0.00034 cm s" 1 


7 




Voltage-gated K + channel sarcolemmal permeability 


0.0004 cm s" 1 


7 


-fk(v,i) 


Voltage-gated K + channel t-system permeability 


0.0001 cm s" 1 


7 



References: 1, Pedersen et al., 2011; 2, Wallinga et al., 1999; 3, model derived to maintain [Na + ]j in physiological range; 4, total conductance taken from 
Pedersen et al. (2011); 5, Pq/^k ratio (Dulhunty, 1979; Pedersen et al., 2005); 6, P^/Pr ratio model derived to obtain membrane potentials in agreement 
with control experiments; 7, Henneberg and Roberge, 1997. Hodgkin-Huxley parameters for voltage-gated Na + and K + channels are from Henneberg and 
Roberge (1997). Note that all intracellular ion concentrations and compartmental volumes in the model are variables, the values of which are determined 
solely by the choice of parameters, as demonstrated in Fig. 2. Extracellular Na + and K + concentrations matched those of the experimental studies, and 
extracellular [Cl~] was chosen to make the solution electroneutral. 
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magnitude and opposite direction, thereby approximat- 
ing stimulation with a KCl-filled microelectrode. 

Calibration to rat skeletal muscle 

The physical conductance parameters summarized in 
Table I were calculated from experimental measure- 
ments of subthreshold sine wave conduction using the 
analytic model of Pedersen et al. (2011). Resting mem- 
brane permeabilities for Na + , K + , and CP were calculated 
from membrane conductance measurements, literature 
values for the ratio Pq/Pk, and model fitting of E m to the 
present experimental data for the P^J Pa ratio. Unless 
otherwise stated, the ionic permeabilities per unit area 
of the t-system membrane and the sarcolemma were as- 
sumed to be identical. The Na + /K + pump was simulated 
using the model of Hernandez et al. (1989), and its 
membrane density was chosen to produce reasonable 
resting values of [Na + ]i. Finally, voltage-gated channel 
permeabilities, distribution, and kinetics were as de- 
scribed by Henneberg and Roberge (1997). 

Model initiation and validation 

A critical advantage of CD modeling is that E^, E t , and 
all intracellular and intra-t-system ion concentrations 
are dependent variables that have unique steady-state 
values for any given set of physical, permeability, and 
pump parameters. Therefore, the model may be initi- 
ated with arbitrary values of the principal modeled vari- 
ables ([Na + ]i, [K + ]i, [Cr]i, [X| [Na + ] t , [K + ] t , [Cr] t , and 
[2] t ). In the present study, the model was run to steady 
state from two markedly different semi-arbitrary sets of 
these variables using a single physiological parameter 
set as listed in Table I. As shown in Fig. 2, the eventual 
steady-state values of all modeled variables were identi- 
cal in each case. This demonstrates that the new model 
fully obeys conservation principles. Furthermore, the 
eventual steady-state values of all variables are physiolog- 
ically reasonable and in agreement with published data, 
thereby providing some validation of the model. This 
procedure was also used to obtain the steady-state values 
in the multiple-shell t-system model (unpublished data) , 
similarly demonstrating a unique steady-state indepen- 
dent of the starting values of the variables. 

MATERIALS AND METHODS 

Animal handling 

Experimental recordings of trains of APs were obtained from sur- 
face muscle fibers of intact EDL muscles from 12-wk-old female 
Wistar rats of our own breed (^230 g). Animals were kept at 21 °C 
and fed ad libitum living under 12:12-h light/ dark conditions. All 
handling and killing of animals followed Danish welfare regulations. 
After dissection, the muscles were incubated in standard Krebs- 
Ringer bicarbonate solution (pH 7.4 at 30°C) containing (in mM) : 
122 NaCl, 25 NaHC0 3 , 2.8 KC1, 1.2 KH 2 P0 4 , 1.2 MgScC 1.3 CaCl 2 , 
and 5.0 D-glucose. All solutions were equilibrated with a mixture 
of 95% 0 2 and 5% C0 2 . 




1000 1500 2000 
Time (s) 

Figure 2. The model cell reaches a single physiologically reason- 
able steady state that is independent of the initial values of any 
modeled variable. All variables are plotted on the left-hand scales, 
except for the potential across the t-system access resistance (£a> 
broken line) , which is plotted on the right-hand scale. The model 
was initiated using physiological parameters (Table I), but with 
two sets of intracellular ion concentrations that were chosen to be 
markedly different from each other and clearly unphysiological, 
although fulfilling a constraint of bulk electroneutrality in each 
compartment. See Table II below for a summary of the initial and 
steady-state values of the major modeled variables. 



Recordings of trains of APs 

Trains of APs were recorded using a two-microelectrode setup de- 
scribed in detail previously (Pedersen et al., 2005). In brief, two 
sharp electrodes (~20 Mfl) were inserted into the same muscle 
fiber, and the train of APs was initiated by injecting an ~3.5-s train 
of constant current pulses (3 ms, 300 nA) into the fiber through 



TABLE II 





From 


legend to Fig. 2 






Variable (units) 


Initial value 1 


Initial value 2 Final value 1 


Final value 2 


[Na + ]; (mM) 


200 


50 


5.84 


5.84 


(mM) 


0 


250 


180.3 


180.3 


[Ct]i (mM) 


0 


50 


7.37 


7.37 


[X]i (mM) 


121.4 


151.7 


108.5 


108.5 


[Na*] t (mM) 


0 


0.01 


157.6 


157.6 


[iT], (mM) 


0.01 


1 


6.38 


6.38 


[Cl] t (mM) 


0 


1.009 


138.0 


138.0 


[Z], (i 1 ) 


-0.01 


-0.001 


-25.4 


-25.4 



Note that the final value of each modeled variable, shown at time <» in Fig. 2, 
is both physiologically reasonable and independent of its initial value. 
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one electrode while the other electrode recorded the membrane 
voltage. To allow the electrodes to remain inserted during AP 
trains, contractile activity of the muscle fibers was reduced by 
including 50 uM of the specific MHC II inhibitor, iV-benzyl-^- 
toluene sulphonamide (BTS) , to the Ringer's solution as described 
by Macdonald et al. (2005) . 



RESULTS 

Our companion study (Pedersen et al., 2011) demon- 
strated that linear electrical circuit models must incor- 
porate a significant t-system luminal resistance to 
accurately model the subthreshold electrical properties 
of rat EDL muscle fibers. It was possible to evaluate the 
physiological significance of such resistances for sarco- 
lemmal AP propagation velocity and t-system excitation. 
This led to an evaluation of the role of Gm regulation 
during repetitive firing of short AP trains (Pedersen 
et al., 2009a). The present study proceeded to a quanti- 
tative analysis of important nonlinear responses not 
amenable to the preceding analytic methods of Pedersen 
et al. (2011) using a new nonlinear model based on 
CD modeling principles (Fraser and Huang, 2004; Fraser 
et al., 2005a,b). First, the frequency dependency of the 
conduction velocity of passively conducted sine waves 
was quantified in the new CD model and compared with 
the previous findings based on linear current-voltage 
relationships (Pedersen et al., 2011). Second, the influ- 
ence of voltage- and time-dependent ion channels within 
both the surface and t-system membranes was clarified. 
Third, this permitted an analysis of t-system ion concen- 
tration changes and their influence upon surface and 
t-system membrane potential and excitability during re- 
petitive AP firing. Experimental microelectrode record- 
ings of membrane potential were compared with the 
results of the CD modeling. 

Passive conduction properties of skeletal muscle fibers 

Fig. 3 A demonstrates the relationship between frequency 
and passive conduction velocity in simulations without 
voltage-gated ion channels under four different degrees of 
separation of the t-system from the true ECF Fiber proper- 
ties were simulated assuming: (a) no R\ (simple cable 
model); (b) a physiological ,RAof72.3 Ocm 2 (homogenous 
t-system model; equivalent to the lumped model of Peder- 
sen et al., 2011); (c) physiological R\ and additional 
smaller series resistances (3.7 mS cm -1 ; Wallinga et al., 
1999) between each of 20 concentric t-system shells (20- 
shell t-system model, equivalent to the distributed cable 
model in Pedersen et al. , 2011); and (d) extremely high R A 
(723 O cm 2 ) simulating near-complete isolation of 
the t-tubules from the ECF (detubulation model) . 

In full agreement with the findings of our companion 
analysis (Pedersen et al., 2011), the homogenous t-system 
model showed faster conduction than the simple cable 
model, particularly at higher frequencies. This was de- 
spite the fundamentally different assumptions and 



approaches underlying the two models. The 20-shell 
t-system model (Wallinga et al., 1999) similarly showed 
faster conduction velocity <1,000 Hz, but slightiy slowed 
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Figure 3. The influence of t-system access and luminal resis- 
tances upon the conduction velocity of a passively propagated 
sine wave and associated t-system potential changes. (A) The fre- 
quency-velocity relationships for model muscle fibers with: no 
resistance between the ECF and the t-system lumen (squares); 
an extremely high access resistance, simulating detubulation (tri- 
angles); a physiological t-system access resistance between a single 
homogenous t-system compartment and the ECF (diamonds); 
and a physiological t-system access resistance and smaller physi- 
ological resistances between each of 20 concentric t-system shells 
(crosses) . Sine waves were driven by a sinusoidal ±300-nA current 
carried by K + and CP ions. (B and C) The underlying changes in 
Em (thick lines), E A (dashed lines), and £ t (o_i9) (lighter lines, from 
Eg, closest to E^, to E lg , furthest from Em) in the 20-shell t-system 
model during stimulation at 50 Hz (B) and 1,000 Hz (C). Note 
that the potential across £i ( o-i9) closely tracks Em at 50 Hz, but lags 
Em in both phase and amplitude at 1,000 Hz. 
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conduction velocity >1,000 Hz, relative to the homoge- 
nous t-system model. It is interesting to note that both 
physiological Ra and intra-t-system series resistances ap- 
pear to enhance the conduction velocity at frequencies 
corresponding to those of the AP upstroke (Pedersen 
et al., 201 1) . Further increases in R A to yield a detubula- 
tion model produced further accelerated conduction at 
lower frequencies (<600 Hz) but slowed conduction at 
all higher frequencies relative to fibers with a physiolog- 
ical R\, regardless of whether additional intra-t-system 
resistances were present. 

Fig. 3 (B and C) demonstrates the changes in t-system 
transmembrane voltage, Et, and the voltage across 
the sarcolemma, Em, during sine wave propagation. In 
full agreement with the analytical solutions of our com- 
panion study, the CD model demonstrates that the volt- 
age division of the intracellular potential between the 
t-system membrane and Ra is highly frequency depen- 
dent. Thus, at low frequency (50 Hz; Fig. 3 B) , E t tracks 
E m relatively closely, indicating that the t-system mem- 
brane then dominates the impedance path from the 
cytosol through the t-system to the ECF At higher fre- 
quency (1,000 Hz; Fig. 3 C), the amplitude of E t is sig- 
nificantly reduced relative to that of Em, corresponding 
to Ra rather than the t-system membrane progressively 
dominating this impedance path. The larger voltage 
drop across R A when compared with that across the 
t-system membrane at high frequencies has two conse- 
quences. First, the t-system membrane is poorly charged 
by high-frequency components of the circuit currents 
during AP propagation. Second, because current and 
voltage do not lag on a resistive component, Ra reduces 
the phase lag between current and voltage at such fre- 
quencies. This reduces the time required for high- 
frequency currents important for AP propagation to 
charge the sarcolemma. This latter point explains 
the steeper frequency-velocity relationship observed 
in models possessing a t-system luminal resistance 
(Fig. 3 A) and furthermore, as in the previous linear 
analysis, suggests that Ra causes the sarcolemmal AP to 
propagate faster. 

The similarities in frequency-velocity relationships 
shown by the CD and linear circuit models (Pedersen 
et al., 2011) occurred despite the distincdy nonlinear 
current-voltage relationships for each ionic current 
and the resultant total currents in the former. These 
nonlinearities result from the differing intracellular 
and extracellular concentrations of the major trans- 
membrane charge carriers, K + and CP, resulting in 
equivalent membrane resistances varying from 2.0 to 
5.3 kfl in the surface membrane and 0.6 to 1.3 kfi in the 
t-system during the 50-Hz sine wave, in each case with 
the lower resistances occurring with greater depolar- 
ization. Thus, the nonlinearity of background cur- 
rent-voltage relationships does not influence passive 
conduction velocity. 



Conduction velocity of actively propagated signals 
in skeletal muscle fibers 

The findings above thus demonstrate that CD model- 
ing successfully reproduces experimental measurements 
and linear modeling of passive frequency-velocity rela- 
tionships in EDL muscle fibers (Pedersen et al., 2011). 
Thus, the simulations of passive conduction demon- 
strate that increased signal frequency, and increased Ra 
at frequencies below ^800 Hz, each produce a greater 
conduction velocity, but also decreased t-system excita- 
tion. However, AP conduction velocities in skeletal muscle 
under physiological conditions depend upon a combina- 
tion of passive signal propagation and active signal regen- 
eration. Furthermore, coordinated excitation-contraction 
coupling along the entire length of a muscle fiber re- 
quires both rapid AP conduction and t-system excita- 
tion. Simulations of regenerative AP conduction were 
therefore performed to: (a) determine the influence of 
active regeneration upon the conduction velocity of 
surface APs; (b) investigate the influence of Ra and the 
presence or otherwise of voltage-gated channels within 
the t-system upon AP conduction velocity and t-system 
excitation; and (c) investigate the influence of the Gm 
changes observed during repetitive AP firing on surface 
and t-system excitability. 

The influence of active signal regeneration upon conduction 
velocity. Fig. 4 demonstrates the importance of active 
signal regeneration to sarcolemmal conduction velocity. 
It displays results from two model fibers with physiologi- 
cal densities of voltage-gated Na + and K + channels after 
stimulation in the central of 99 compartments into 
which the fiber long axis was divided. Fig. 4 A compares 
voltage responses in six adjacent compartments 3-3.5 mm 
from the initial stimulus along each model fiber. The 
solid lines depict results from a model fiber with a physio- 
logical density of voltage-gated ion channels. The bro- 
ken lines depict the passive conduction remaining when 
the permeabilities of the voltage-gated ion channels in 
all compartments beyond 3 mm from the stimulus were 
fixed at their resting values. The latter shows that signal 
amplitude steeply declines with distance, reflecting an 
absence of active signal regeneration. Furthermore, the 
increased spacing of the peaks with distance in the pas- 
sive conduction model (marked by vertical broken 
lines) indicates a marked deceleration of the passively 
propagated signal with distance from its source, relative 
to the actively propagated signal. 

Fig. 4 (B-D) goes on to examine the detailed effects 
of Pns upon the conduction characteristics. Fig. 4 B (inset) 
shows AP traces resulting from setting the maximum 
permeability of the voltage-gated Na + channel (-PNa(max)) 
at 200% (i), 100% {it), 50% {Hi), or 30% (iv) of the 
value expected under control conditions (Henneberg 
and Roberge, 1997). The AP magnitude increased and 
the time to peak decreased with increasing PNa(max)- 
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The resulting sarcolemmal conduction velocity of ac- 
tively propagated APs (Fig. 4 B, squares) then showed a 
log-linear dependence upon -PNa(max), similar to the rela- 
tionship between passive conduction velocity and sine 
wave frequency (Fig. 3). This contrasted with the shal- 
lower dependence on PNa(max) of the velocity of passively 
conducted AP waveforms over the first 100 pm ahead of 
the active AP wave front. This is consistent with the 
highest frequency, fastest conducting components of 
the AP having the shortest length constants (Pedersen 
etal., 2011). 

Fig. 4 C summarizes a power spectral analysis of the 
AP traces in the inset of Fig. 4 B obtained by their 



Fourier transformation to permit a comparison of 
their frequency content and the conduction velocity 
of the corresponding APs. It clearly correlated an in- 
creased PNa(max) with a larger content of high-frequency 
components in the resulting AP waveforms. This could 
be compared direcdy with the frequency dependence 
shown by the propagation velocity of passive sub- 
threshold sinusoidal voltage changes, through identi- 
fication of an equivalent frequency for the AP equal 
to the sine wave frequency, whose conduction velocity 
matched that of the AP. As passive frequency compo- 
nents with conduction velocities slower than that of 
the AP cannot contribute to the AP propagation velocity, 




o Relative P Na( max ) 

Frequency (Hz) ^ 1 1 J 1 1 ± Q 

Figure 4. Comparison of passive conduction with active AP regeneration. The model with homogenous t-system was stimulated in the 
central of the 99 compartments at 2 s. (A) £„, over time for segments 30-35 (3-3.5 mm from the stimulus) from two simulations: in the 
first simulation (solid lines), active propagation persisted throughout all fiber segments; in the second (broken lines), conduction block 
was modeled in compartments 31 and beyond by fixing the permeability of the voltage-gated ion channels to their resting values. This 
allowed a direct comparison of active and passive signal propagation. The vertical broken lines mark the peaks of these passively con- 
ducted AP waveforms. It can be seen that the peaks of the passively conducted signals eventually lag the peaks of the actively propagated 
APs. (B) The measured peak-to-peak conduction velocity in similar simulations of conduction block, over a range of values of P Na(max) , 
for actively propagated signals (squares); and for the first 100 pm of passive conduction from the final point of active regeneration (tri- 
angles). (Inset) Four traces of APs with different relative values of iNA(max) 0.3; », 0.5; in, 1; and iv, 2). The frequency content of these 
traces of APs was analyzed using Fourier analysis, and their power spectra are shown in C. Given that APs propagate via the spread of 
local circuit currents in front of the AP waveform, some elements of such passive current flow must be faster than or as fast as the active 
propagation. Thus, from the relation of frequency and passive conduction velocity of sinusoidal currents displayed in A, the frequency 
of a sinusoidal current required to match the AP propagation velocity was determined for the range of P Na(max) values investigated. These 
frequencies were referred to as equivalent frequencies. (D) This illustrates that the equivalent frequency and the active AP conduction 
velocity were linearly related and, furthermore, that the equivalent frequency and iNa(max) were log-linearly related. 
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this equivalent frequency therefore identifies the 
lowest frequency component of the circuit currents 
that could possibly contribute to the AP propagation 
velocity. Fig. 3 A was used to extract the equivalent fre- 
quencies for the different conduction velocities ob- 
tained at the different PNa(max) values. Fig. 4 D plots the 
resulting equivalent frequencies against either AP 



conduction velocity or -PNa(max)- The relationships 
were linear for the conduction velocity and log-linear 
for i°Na(max) ranging from 0.3 to 10 times the control 
value. Collectively, these findings clearly illustrate the 
influence of -PNa(max) upon the high-frequency compo- 
nents of the AP, and through the latter upon AP con- 
duction velocity. 
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Figure 5. The influence of t-system excitability upon the t-system response to a sarcolemmal AP. The response of i4(„) (thin lines) and E A 
(broken lines) to changes (thick lines) during an AP is depicted for a model fiber with a t-system simulated as a single homogenous 
compartment (A and B) and for a model fiber with a 20-shell t-system (C and D). In A and C, the t-system lacks voltage-gated ion chan- 
nels and therefore the voltage response is entirely passive. B and D demonstrate the response of a fiber with physiological voltage-gated 
ion channel densities within the t-system. Surface £„, changes are similar in each case, but the resulting tubular excitation is enhanced 
by the presence of voltage-gated ion channels within the t-system. (E) A summary of peak changes in fiber core during the passage 
of a sarcolemmal AP. (F) Surface conduction velocity for fibers with (squares) and without (triangles) voltage-gated ion channels in the 
t-system, in the homogenous t-system model (dashed lines) and in the 20-shell t-system model, for model fibers with a range of values of 
R A (solid lines). The vertical line in each case denotes the value of Ba used for the simulations in A-D. Increases in R\ reduce the passive 
electrical response of the t-system to the sarcolemmal AP (E) while accelerating sarcolemmal conduction velocity (F). A physiological 
value of R A thus necessitates voltage-gated channels to produce a full t-system voltage response (E). Such channels have no direct effect 
on sarcolemmal conduction velocity (F) . 



Fraseretal. 105 



Control 



Phase I 



Phase I 



I I I I I 



JJJJ 
I I I I I I 



.Mil 



II 



B 




-60 -, 


> 




£ 


-65 - 






.tent 


-70 - 


0. 




a> 


-75 - 


to 




-Q 

E 






-80 - 




-85 - 



D 




> Control 
i Phase I 
i Phase II 



I 1 

1 11 1 12 

Time (s) 




I 

1.13 



"I 1 1 1 

1.200 1.202 1.204 1.206 

Time (s) 




— I 1 1 1 

2100 2 102 2 104 2 106 

Time (s) 



Control 



Phase I 





Phase I 



— Sarcolemma 

— t-system shells 




Time (ms) 



Time (ms) 



Control 




Time (ms) 



Phase I 



— Sarcolemma 

— t-system shells 




Phase I 



— Sarcolemma 

— t-system shells 




1 

4 

Time (ms) 



Figure 6. Effects of G M regulation during Phase I and Phase II on simulations of excitation of sarcolemmal APs, their propagation, and 
t-system excitation. (A) Current flow during endplate excitation of sarcolemmal APs was mimicked by simulating injections of 1-ms square KC1 
currents in the central segment of the model. To determine the effect of Cm regulation on the current required to excite an AP, a range of 
different current amplitudes was applied. (A) The membrane potential of the central segment of the model where excitation occurred (the 
current injections are presented below). Membrane permeabilities of the model for K + and Cl~ were adjusted to mimic G M during control 
conditions (left), during Phase I (middle), and during Phase II (right). (B) The membrane potential response to a current injection that 
did not elicit an AP under any of the conditions. Clearly, Gm regulation markedly affects the recovery of membrane potential after cessation 
of current injection. (C) Membrane potential under the three conditions of control, Phase I, and Phase II during the injection of a current 
that only triggered an AP during Phase I. (D) The membrane potential when the current amplitude was just sufficient to elicit an AP in 
Phase II. (E) Sarcolemmal AP (black line) and the membrane potential in the 20 t-system shells 3 mm from the point of AP excitation dur- 
ing control conditions (left), Phase I (middle), and Phase II (right). (F) Similar to E, but with a 50% reduction in maximum permeability 
for voltage-gated Na + channels. 
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The effect of t-system excitability upon surface conduction 
velocity and t-system excitation. The simulations depicted 
in Fig. 4 demonstrate that, over anything but the short- 
est distances, active signal regeneration by sarcolemmal 
voltage-gated ion channels enhances conduction veloc- 
ity. The simulations that follow now analyze the effects 
of the voltage-gated ion channels additionally present 
in the t-system membrane upon t-system excitation and 
surface conduction velocity. 

Fig. 5 compares passive (A and C) and active (B and 
D) t-system responses to a surface AP in a model with a 
single homogenous t-system compartment (A and B) 
and in a model with a t-system subdivided into 20 
concentric shells of equal thickness (C and D). In each 
case, changes in E^ lag changes in E^. In the single- 
compartment t-system model, the magnitude of the 
change in E t ( A£ t ) is increased from +73.5 mV in the ab- 
sence of voltage-gated ion channels within the t-system 
(Fig. 5 A) to +109.8 mV with active regeneration (B). 
However, findings from the 20-shell t-system model sug- 
gest that the homogenous t-system model underestimates 
the importance of active regeneration within the t-system. 
Thus, the 20-shell model shows that the E t change in the 
center of the fiber (AE^ig^) is increased from +40.9 mV 
in the absence of voltage-gated ion channels within the 
t-system (Fig. 5 C) to +118.0 mVwith active regeneration 
(D) . This contrast between the two t-system models is of 
some relevance, as A£i in the passive single-compartment 
t-system model would be sufficient for activation of 
excitation-contraction coupling, whereas AE^igj in the 
passive 20-shell t-system model is well below this activa- 
tion threshold (Gallant and Jordan, 1996). 

Fig. 5 (E and F) demonstrates the influence of Ra 
magnitude upon the peak E t change during an AP (E) , 
and upon the sarcolemmal AP conduction velocity (F) . 
In the homogeneous model (Fig. 5, E and F, broken 
lines) , an increase in R\ decreases the passive change in 
Ei during AP in fibers without voltage-gated ion chan- 
nels in the t-system (E, triangles), but accelerates the 
sarcolemmal AP conduction velocity (F, triangles). Sim- 
ilarly, the presence of a series resistance separating suc- 
cessive t-system shells in the 20-shell t-system model 
(Fig. 5, E and F, solid lines) produces a further modest 
acceleration of the sarcolemmal AP (F), while greatly 
decreasing the passive t-system voltage response (E), 
relative to that of the homogenous t-system model. How- 
ever, in the presence of voltage-gated ion channels in 
the t-system (Fig. 5, E and F, squares), the change in E^ 
during an AP is >100 mV in both the homogenous and 
20-shell t-system models, even for values of Ra up to 
three times the physiological value. 

Thus, these findings demonstrate that t-system lumi- 
nal resistance enhances surface AP conduction at the 
expense of a reduced passive excitation of the t-system. 
However, even in the presence of a significant access re- 
sistance, voltage-gated ion channels within the t-system 



permit full t-system excitation during an AP. Finally, 
voltage-gated channels in the t-system do not contribute 
to the velocity of surface AP conduction. 

The influence ofG M regulation in active muscle on sarcolem- 
mal and t-system excitability. Despite the advantages of 
the analytical solution in generating a detailed mecha- 
nistic insight into the relationship between Gm and ex- 
citability (Pedersen et al., 2011), the approach was 
limited to linear analysis. The present study therefore 
extended the work to determine the influence of Gm 
changes upon muscle excitability in the presence of 
voltage-gated ion channels in sarcolemma and in the 
t-system. The simulations test the prediction that AP ex- 
citation at the neuromuscular junction and the t-system, 
but not along the sarcolemma, should be particularly 
sensitive to alterations in Gm, owing to sarcolemmal AP 
propagation being conveyed by high-frequency ele- 
ments of the circuit currents for which the impedance 
of a sink membrane region would be relatively insensi- 
tive to Gm- 

Thus, Fig. 6 shows the influence of Phase I and 
Phase II Gm changes (see Fig. 9 C) on the excitation of 
the sarcolemma AP. During Phase I, it was possible to 
trigger an AP with a current that was slightly lower than 
the current that triggered an AP under control condi- 
tions, whereas during Phase II, a substantially larger cur- 
rent was required. Thus, sarcolemmal AP excitation is 
facilitated by Phase I Gm regulation and markedly im- 
paired by Phase II G M regulation, even when excitation 
proceeds with currents of short durations comparable 
to the endplate current. 

Fig. 6 B shows the subthreshold responses under each 
condition. Gm strongly influences the downstroke of the 
membrane potential response, exerting smaller influ- 
ences on its upstroke and magnitude. The importance 
of this influence on the downstroke is demonstrated in 
Fig. 6 C, in which it is shown that for near-threshold 
stimuli, the AP upstroke occurs after the cessation of 
the 1-ms stimulus. Fig. 6 D additionally demonstrates 
that APs excited during Phase II were characterized by 
a reduced amplitude and slower upstroke as compared 
with APs during Phase I and under control conditions. 

Next, sarcolemmal AP propagation and t-system exci- 
tation were simulated in the model during control con- 
ditions, Phase I, and Phase II. Fig. 6 E shows the 
sarcolemmal AP (black lines) and the membrane volt- 
age in the t-system shells (red, blue, and green lines) at 
3 mm from the point of excitation. Under control con- 
ditions, the simulated APs propagated at 2.03 m s~ , 
whereas during Phase I and Phase II, it propagated at 
2.09 and 1.67 m s , respectively. Thus, Phase I caused 
only a 3% increase in propagation velocity, whereas 
Phase II reduced the propagation velocity by 18%. 
Under conditions when the maximum permeability of the 
voltage-gated Na + channels was reduced to 50%, the 
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sarcolemma propagation velocity dropped to 1.52 m s~ 
under control conditions (Fig. 6 F) . The increase in sar- 
colemmal AP propagation velocity during Phase I be- 
came marginally more pronounced (5%), whereas dur- 
ing Phase II, the sarcolemmal AP propagation velocity 
dropped by 36%. There were also clear-cut effects of Gm 
regulation on t-system excitation and t-system AP wave- 
form. Thus, the amplitude of the AP in the innermost 
t-system shell became 12 and 14% larger during Phase I, 
whereas during Phase II, it was reduced by 19%. With 
reduced -PNa(max), the t-system shells did not actually pro- 
duce an AP, despite the existence of fully propagating 
sarcolemmal AP during Phase II. The latter finding sug- 
gests that Gm regulation during Phase II can cause loss 
of t-system excitation when voltage-gated Na + channels 
have been partly inactivated, as can occur through slow 
inactivation of these channels during repetitive AP fir- 
ing. This suggests that Gm regulation could contribute 
to muscle fatigue through reductions in both endplate 
and t-system excitability. 

Changes in t-system ionic homeostasis and membrane 
potential during AP trains 

Skeletal muscle function in vivo generally involves se- 
quential as opposed to lone AP firing. The approach 
used in this analysis permitted study of ion concentra- 
tion and membrane potential homeostasis associated 
with such repetitive firing. It is generally established 
that repetitive AP firing in skeletal muscle fibers mark- 
edly elevates t-system luminal [K + ] (Aimers, 1980; 
Wallinga et al., 1999; Shorten and Soboleva, 2007). It 
has also been speculated that these could depolarize 
the t-system membrane sufficiently to produce signifi- 
cant Na + channel refractoriness and thereby reduce 
t-system excitability. Thus, tubular K + accumulation 
has been implicated in muscular fatigue (Sj0gaard, 
1990). Experimental studies report such accompany- 
ing depolarizations of the resting membrane potential 
(Freygang et al., 1964). 

However, experimental explorations of the possible 
tubular K + changes, depolarization phenomena, the 
relationship between them, and the extent to which 
the voltage changes might provide a means for assess- 
ing tubular K + homeostasis during AP firing are pre- 
cluded by real-time assessments for [K + ] t not being 
available. Nevertheless, CD modeling provided an ap- 
proach to such an analysis, as it calculates the concen- 
trations of all ions in all compartments. The present 
studies determined the extent to which depolarization 
of the resting membrane potential during an AP train 
in skeletal muscle fibers might result from t-system K + 
accumulation, through comparing experimental re- 
cordings with simulations of AP trains using the 20- 
shell t-system model. 

Fig. 7 A shows representative experimental recordings 
of AP trains with AP-firing frequencies of 6, 15, and 30 Hz 



in EDL muscle fibers. In each case, the resting membrane 
potential developed a depolarization over the first M).5 s 
of AP firing that then remained relatively stable despite 
continuation of the AP trains for 3.5 s. The magnitude of 
this depolarization increased with the AP-firing frequency. 
Thus, on average, the maximum baseline depolarization 
was 0.97 ± 0.15 mV at 6Hz(b = 13), 2.08 ± 0.12 mV at 
15 Hz (n = 23), and 4.4 ± 0.33 mV at 30 Hz (n = 19). On 
cessation of AP firing, the resting membrane potentials re- 
covered over — 1 s to their prestimulation value. Note that 
the resting membrane potential during AP firing was re- 
corded immediately before each AP was triggered. 

Fig. 7 B shows sarcolemmal potentials during simu- 
lated AP trains of the same duration and frequency as in 
the experimental records depicted in Fig. 7 A. Fig. 7 C 
then compares the experimental and model records of 
sarcolemmal depolarization, measured as the difference 
between the potential immediately before each stimula- 
tion and that immediately before the first AP. The CD 
model successfully replicated the magnitudes of the 
resting membrane potential depolarizations seen dur- 
ing the experimentally recorded AP trains. 

To determine whether the depolarization of the rest- 
ing membrane potential was associated with alterations 
in t-system [K + ] , Fig. 7 D shows t-system [K + ] for each of 
the 20 t-system shells during the AP trains. AP firing led 
to marked elevations in [K + ] t , particularly in the deeper 
regions of the t-system. At all frequencies, [K + ] t reached 
a steady level after — T s of AP firing and rapidly recov- 
ered after AP firing. The magnitude of change in [K + ] t 
increased with frequency. Thus, mean t-system [K + ] t in- 
creased from 5.3 mM before AP firing to —6.8 mM at 
6 Hz, -8.8 mM at 15 Hz, and -41. 2 mM at 30 Hz. 

The cause of t-system [K + ] accumulation was then fur- 
ther explored. First, simulation of fibers without t-system 
voltage-gated Na + channels showed a depolarization of 
only 0.3 mV when stimulated at 15 Hz for 4 s, in compari- 
son with the 2.1-mV depolarization in modeled fibers 
with normal t-system ion channel densities. This re- 
flected a much reduced accumulation of K + within the 
t-system (increase of [K + ] t of around 0.5 mM compared 
with 5.9 mM under normal conditions). Similarly, there 
was an absence of resting potential depolarization in a 
model of detubulation (unpublished data). Therefore, 
the model suggests that t-system AP firing is responsible 
for the K + accumulation and baseline membrane poten- 
tial depolarization seen during experimental AP trains, 
and therefore might be used experimentally to assess 
whether sarcolemmal APs trigger APs in the t-system. 

The simulations demonstrated much smaller sarco- 
lemmal and tubular membrane depolarization than 
would be expected simply from the alterations in the 
equilibrium potential of K + that the accumulation of K + 
in the t-system would induce. This suggested a role for 
CP in stabilization of the resting potential during re- 
petitive AP firing. 
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The role of Cr conductance on t-system K + homeostasis 
and the membrane potential during AP trains 

Fig. 8 explores the possible role of CP in simulations 
and experimental measurements of repetitive AP fir- 
ing using solutions with reduced CP concentration. 



The replacement of extracellular CP with a monovalent 
membrane-impermeant anion produced a marked in- 
crease in the magnitude of the depolarization observed 
during repetitive AP firing in both the experimental re- 
cordings (Fig. 8, A and B) and the simulations (Fig. 8 C). 
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Figure 7. Comparison of AP trains recorded experimentally with those simulated using the CD model. Trains of APs show a similar 
baseline depolarization of the membrane potential in experimental recordings from rat EDL muscle fibers (A) and in model simula- 
tions (B) when stimulated at 6, 15, or 30 Hz. The resting membrane potential during AP trains was evaluated as the membrane potential 
immediately before injection of excitatory current through the current electrode. At each frequency, the depolarization of resting mem- 
brane potential (C) was similar in the experimental (open circles; mean ± SEM) and model records (solid lines). In the model (D), this 
depolarization resulted from K + buildup in the t-system that was greatest in the innermost t-system shell (uppermost traces) , and least in 
the outermost shell (lowest traces) . 
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In previous studies (Pedersen et al., 2009a, b), Gq was 
measured in EDL muscle fibers at the different CP 
concentrations also used in the present study. 

Under control conditions Gci and Gi^were 1,314 and 
144 pS/cm 2 , respectively, giving a Gk/Gq of 0.11. Re- 
ducing extracellular CP to 80 mM increased the depo- 
larization magnitude by around 50% and increased 
Gr/Gci to 0.28. Further reduction in CT to 50 mM in- 
creased Gg/Ga to 0.66 and the depolarization by 143%. 
Entirely withdrawing CP increased the depolarization 
magnitude by 543%, despite the necessary addition of a 
small dose of TTX (10~ 8 M) in the solutions to avoid 
myotonic activity. 

These experiments suggest a highly nonlinear rela- 
tionship of the depolarization during AP firing to 
Gk/ Gci, in which moderate reductions in Gci, similar to 
those during Phase I (Gk/ Gci of around 0.5) , only mod- 
erately increase the depolarizations, whereas further 
Gci reductions markedly increase the depolarizations dur- 
ing AP firing, potentially giving rise to myotonic activity. 

Simulations of AP trains with different CP concentra- 
tions gave similar findings. Thus, Fig. 8 C shows that for 
control conditions and for 80 and 50 mM CP, the mag- 
nitude of depolarization during the simulated trains 
closely agreed with the experimental observations. The 
model further showed that the increased depolarization 



with reduced CP was associated with minor increases 
in concentrations of t-system K + during the AP trains 
(Fig. 8 D) . In accordance with the experiments in CP-free 
conditions where TTX had to be included to prevent 
myotonic AP firing, the simulations in the absence of 
CP resulted in myotonic activity. It was therefore neces- 
sary to reduce the maximum permeability of voltage- 
gated Na + channels in the model to perform simulations 
without CP. Fig. 8 C shows that the depolarization dur- 
ing the simulated AP train in CP-free conditions caused 
a marked rise in the depolarization that, however, was 
not as marked as in the experimental recordings. In 
CP-free conditions, the rise in t-system K + during the AP 
trains was similar to the rise at 50 mM CP. 

Gm regulation during repetitive AP firing has been 
shown to produce first a decrease in CP conductance 
(Phase I), and then a large increase in K + and CP con- 
ductances (Phase II) (Pedersen et al., 2009a, b). The in- 
fluence of these changes on the magnitude of baseline 
membrane potential depolarization and t-system ionic 
homeostasis during repetitive AP firing is demonstrated 
in Fig. 9. Fig. 9 A shows the first, 40th, and 80th trains of 
49 APs at 15 Hz and, on a larger scale, the voltage response 
to a square test current pulse injected in between the 
AP trains. This demonstrates an increased voltage re- 
sponse immediately before the 40th AP train, illustrating 




t 1 1 1 1 1 i 1 1 1 1 1 i 1 1 1 1 1 

234567 234567 234567 

Time (s) Time (s) Time (s) 



Figure 8. The effect of extracellular CP on the magnitude of the depolarization during AP trains. Two electrodes were inserted into 
the same muscle fiber, and 3.5-s-long 15-Hz trains of AP trigger pulses were injected. Experiments were conducted using extracellular 
solutions without CI" (n=20) orwith different CP concentrations of 127 (control; n= 11), 80 (« = 23),or50 mM (n= 14). (A) Repre- 
sentative recordings from different fibers at the four CP concentrations. (B) The average change in resting membrane potential during 
the AP trains at the four CP concentrations. Also presented is the Gci at the different conditions as reported by Pedersen et al. (2009a) 
and the resulting Gk/ G a ratios. (C) The experimental observations of AP trains at the different CP concentrations were simulated, and 
the change in resting membrane potential was determined. (D) The K + concentration in the deepest t-system compartment during the 
AP trains under the four conditions. 
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the Phase I Gm decrease, and a decreased voltage re- 
sponse before the 80th AP train, illustrating the Phase 
II Gm increase. Fig. 9 B demonstrates that the Phase I Gm 
decrease results in an approximately twofold increase in 
baseline depolarization, similar to the increase in base- 
line depolarization seen in low CP solutions (Fig. 8) . 
Yet, a smaller increase in baseline depolarization com- 
pared with initial conditions was also recorded during 
the Phase II increase in Gm- 

The precise changes in Gm (Fig. 9 C) and G$J G a ra- 
tios (D) were obtained from Pedersen et al. (2009a) to 
calibrate the model. Thus, Fig. 9 E demonstrates simu- 
lations showing a similar pattern of a greatly increased 
baseline depolarization during Phase I and a smaller in- 
crease in baseline depolarization during Phase II Gm 
changes. However, Fig. 9 F demonstrates that although 
both Phase I and Phase II increase the baseline depolar- 
ization, these Gm changes have opposite influences on 
t-system K + accumulation. Thus, t-system K + accumula- 
tion was increased in Phase I by ^6 mM compared with 
4.5 mM in control conditions, whereas it only was in- 
creased by 2 mM in Phase II. In addition, larger excur- 
sions in [K + ] t were seen during each AP in Phase II. Finally, 
Fig. 9 G shows the influence of membrane CP permea- 
bility on baseline depolarization and t-system K + accu- 
mulation. At low values of Pci, K + accumulation and 
baseline depolarization are large, as shown in Fig. 8. Be- 
tween values of Pci of 10 cm s~ and 10~ 6 cm s , there 
is a log-linear relationship between P C i and baseline de- 
polarization, but a smaller influence of Pa on t-system K + 
accumulation. Thus, in this range that includes the Pa 
seen in Phase I, the predominant effect of a decrease in 
Pci is to shift the membrane potential closer to E K . 
However, an increase in P a above normal values, as seen 
in Phase II, produces a sharp reduction in t-system K + 
accumulation and yet an increase in the magnitude of 
the baseline depolarization. This latter effect reflects 
depolarization of £ci because of an increase in intracel- 
lular [CP] as CP currents, rather than K + currents, pro- 
vide the majority of repolarization after each AP. Note 
that further increases in P a prevent t-system excitation. 

DISCUSSION 

The t-system of skeletal muscle fibers both influences 
and is influenced by the passage of a sarcolemmal AP. It 
has a capacitance that is considerably greater than that 
of the sarcolemma itself (Falk and Fatt, 1964; Gage and 
Eisenberg, 1969), and therefore its charging could po- 
tentially slow the surface conduction velocity. It also 
possesses voltage-gated ion channels that underlie the 
t-system AP (Bezanilla etal., 1972; Posterino etal., 2000; 
Pedersen et al., 2004), which is reflected in recordings 
of the surface AP as a prominent after-depolarization 
(Gage and Eisenberg, 1967) . Such t-system regenerative 
electrical activity is thought to result in an accumulation 



of K + in the t-system lumen during repetitive firing of 
APs (Wallinga et al., 1999; Sejersted and Sj0gaard, 2000) . 

It has been proposed that such K + accumulation could 
induce loss of muscle excitability through depolarization 
and consequent inactivation of voltage-gated Na + chan- 
nels, thereby contributing to muscle fatigue (Sejersted 
and Sj0gaard, 2000). Key to this theory on the etiology 
of muscle fatigue are observations of elevated extracel- 
lular K + in plasma (Nordsborg et al., 2008) and in the 
interstitial space of contracting muscle as obtained 
using microdialysis probes (Mohr et al., 2004). However, 
for two reasons such measurements of extracellular K + 
confer very little information on muscle excitability and 
their membrane potential per se. First, the membrane 
potential depends heavily on the actions of other ions, 
their permeabilities, and other transport systems. Sec- 
ond, it has been shown that several mechanisms are ac- 
tivated in the active muscle that can either reduce or 
enhance the force-depressing actions of activity-induced 
elevation in extracellular K + . Thus, muscle activity is asso- 
ciated with elevated Na + ,K + pump activity (Overgaard 
and Nielsen, 2001) and changes in Gm (Pedersen 
et al., 2009a,b) that might also markedly affect the excit- 
ability in working muscle. To deal with the emerging 
complexity of skeletal muscle excitability, it is desirable 
to develop mathematical models that will allow multiple 
parameter changes to be evaluated in combination. 

In this study, a model was developed that builds on a 
realistic t-system geometry as determined using the ana- 
lytic approach in our companion paper (Pedersen et al., 
2011) to evaluate the role of Gm regulation upon excit- 
ability and AP conduction velocity in a system with non- 
linear voltage- and time-dependent conductances, and 
to describe and quantify the interrelationships between 
Gm, t-system ionic homeostasis and the membrane po- 
tential during repetitive AP firing. 

The new multiple-compartment model was based on 
the CD model of Fraser and Huang (2004) and Fraser 
et al. (2005a,b) , previously used to investigate ionic and 
volume homeostasis in amphibian skeletal muscle. The 
model included background conductances, t-system re- 
sistances, and a longitudinal cable structure calibrated 
to rat skeletal muscle using the analytic model of Pedersen 
et al. (201 1) , and voltage- and time-dependent conduc- 
tances to allow simulation of regenerative AP firing 
(Wallinga et al., 1999). The ion concentrations within 
all compartments of the model were shown to possess 
unique history-independent steady-state values deter- 
mined by the physical and ion channel parameters of 
the model. This allowed for the use of the model to pre- 
dict ion concentrations within the t-system both at rest 
and during repetitive activity, as well as sarcolemmal 
and t-system AP conduction. 

The model was first used to simulate the passive con- 
duction properties of rat skeletal muscle in the absence 
of voltage-gated ion channel activity to extend the work 
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of Pedersen et al. (2011) into a model system with physi- 
ological nonlinear background current-voltage rela- 
tionships. This yielded similar results to those obtained 



from the model system with linear current-voltage rela- 
tionships. Thus, the passive conduction velocity of a si- 
nusoidal input current was greater for higher frequency 
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Figure 9. Effects of G M regulation during Phase I and Phase II on the membrane potential depolarization and t-system K + concentration 
during trains of APs. (A) Representative experimental recordings of the membrane potential from a muscle fibers that was repeatedly 
activated to fire 49 APs at 15 Hz every 7 s. Recordings show the first (left) , the 40th (middle) , and the 80th (right) AP train. To highlight 
the regulation of G M during AP firing, the membrane potential responses (underlined above) to the small current injections in between 
the AP trains have been highlighted below. As previously reported (Pedersen et al., 2009a), the onset of AP firing leads to a reduction 
in G M , as shown by an enlarged membrane potential deflection during the current injection (Phase I), whereas prolonged AP firing 
leads to a marked rise in G M (Phase II), as revealed by a reduced membrane potential deflection. To determine the effect of such G M 
regulation in active muscle fiber on the depolarization of the resting membrane potential during the AP trains, the membrane poten- 
tial before AP firing was subtracted from the potential as recorded 55 ms after every AP. (B) The average change in resting membrane 
potential during the AP trains during control conditions (first AP train), during Phase I (40th train), and during Phase II (80th train). 
(C) Gm during AP firing as originally reported by Pedersen et al. (2009a) . (D) Changes in the composite conductances for K + (Gk) and 
Cl~ (Gci) that underlie the changes in Gm in C. (E) The depolarization of the resting membrane potential during the AP trains was 
determined in simulated AP trains in a similar way that this depolarization was assessed in the experimental recordings in B. (F) The K + 
concentration in the deepest t-system compartment during the simulation of AP trains during the three conditions. (G) The magnitude 
of depolarization of the resting membrane potential during simulated 15-Hz trains (circles) and the corresponding t-system K + in the 
deepest t-system shell (triangles) for a range of CP membrane permeabilities. Typical values for CV permeabilities for the three condi- 
tions have been indicated. 
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signals, and this frequency-velocity relationship was 
greatly enhanced by the presence of a significant access 
resistance (R\) between the t-system and the extracellu- 
lar space surrounding the fiber. Furthermore, the pres- 
ence of additional series resistances within the t-system 
was shown to steepen this frequency-velocity relation- 
ship at higher frequencies (above ^300 Hz). 

The study then modeled the influence of active sig- 
nal regeneration by voltage-gated ion channels upon 
surface conduction velocity and t-system excitation. 
This demonstrated that entrance of a propagating AP 
into a region lacking voltage-gated ion channels pro- 
duced a sharp deceleration of conduction velocity 
over distance. Thus, the presence of AP regeneration 
by voltage-gated ion channels allowed actively propa- 
gated APs to show faster peak-to-peak conduction than 
did passively propagated AP waveforms, even over very 
short distances (100 um or less), despite the clear ne- 
cessity for passive currents to travel ahead of the AP 
waveform. Fourier transformation demonstrated that 
high-frequency components of passively conducted 
AP waveforms declined steeply over distance, thereby 
accounting for the reduced conduction velocity. In 
contrast, active regeneration by voltage-gated ion 
channels was shown to maintain the high-frequency 
components and thereby maintained a high conduc- 
tion velocity. Thus, the analysis predicted a log-linear 
relationship between the maximum total Na + permea- 
bility (j°Na(max)) and AP velocity over a range of values 
of iNa(max) from 0.3 to 10 times the normal physiologi- 
cal value. Fourier analysis then showed that this influence 
of ^Na(max) upon conduction velocity could be attributed 
to a greater content of high-frequency components in 
the AP and therefore in the local circuit currents that 
convey AP propagation. 

Fast AP conduction is necessary but not sufficient 
to achieve synchronous contraction along the entire 
length of a muscle fiber. There is also a requirement for 
the surface AP to excite the entire t-system sufficiently 
to activate dihydropyridine receptors and thereby initi- 
ate excitation-contraction coupling (Huang and Peachey, 
1989). The present study therefore explored the in- 
fluence upon t-system excitation of the passage of a 
sarcolemmal AP for several different t-system models. 
It demonstrated that a model without a t-system luminal 
series resistance did not require voltage-gated ion chan- 
nels within the t-system in order for a surface AP to pro- 
duce sufficient t-system depolarization (to ^0 mV) for 
excitation-contraction coupling to occur, despite a 
physiological -Ra- In contrast, a t-system simulated as 20 
concentric shells with small physiological resistances be- 
tween each shell demonstrated a much reduced re- 
sponse (to a peak of only approximately —40 mV in the 
fiber core) to the passage of a surface AP in the absence 
of voltage-gated ion channels within the t-system. How- 
ever, the presence of voltage-gated ion channels within 



the t-system allowed for the propagation of a full-sized 
AP through the t-system, producing a peak response of 
approximately +40 mV in the core of the fiber. The pres- 
ence or otherwise of voltage-gated Na + channels within 
the t-system had no influence on the conduction veloc- 
ity of the sarcolemmal AP. Thus, it appears clear that fast 
AP propagation requires a significant Ra, but a high R A 
reduces the passive t-system excitation produced by a 
surface AP, such that regeneration by voltage-gated ion 
channels within the t-system is required. 

Our companion study (Pedersen et al., 2011) addi- 
tionally demonstrated that Ra has an important influ- 
ence on the relationship between G M and the response 
of the sarcolemma and the t-system membrane to sub- 
threshold signals. Thus, it predicted that fiber excitation 
at the neuromuscular junction and t-system excitation 
are particularly sensitive to regulation of Gm, compared 
with sarcolemmal propagation. This prediction arose 
from the property that sarcolemmal AP propagation is 
conveyed by high-frequency elements in the circuit cur- 
rents, and the membrane impedance at such frequen- 
cies is largely insensitive to Gm- Despite the advantages of 
the analytical solution in generating the detailed mecha- 
nistic insight into the relationship between Gm and excit- 
ability, the approach was limited to a linear analysis of 
subthreshold responses. It was therefore useful to ex- 
tend the predictions generated in our companion paper 
using a model system simulating regenerative AP firing. 

During repetitive AP firing, there is first an M50% re- 
duction in Gm (Phase I) , and then after prolonged activ- 
ity, there is a four- to fivefold increase in Gm (Phase II) , 
in comparison to control conditions (Pedersen et al., 
2009a,b). The simulations demonstrated that, for 1-ms 
square current pulses simulating an endplate potential, 
Phase I produced a minor decrease in the current re- 
quired to excite an AP, whereas Phase II required a near 
doubling of this current. This suggests that the sarcolem- 
mal AP excitation is somewhat facilitated by Phase I and 
markedly reduced during Phase II Gm regulation, even 
though currents of short duration, comparable to the end- 
plate current, must predominantly flow across the mem- 
brane capacitance. Phase II Gm changes were also shown 
to produce a reduction in the upstroke velocity of the sur- 
face AP sufficient to reduce conduction velocity by ^18%, 
and a reduction in the magnitude of the t-system AP. 
Additional reduction of PNa(max) — simulating activity in- 
duced slow inactivation of the voltage-gated Na + channel — 
produced a failure of t-system AP generation during 
Phase II Gm changes, despite persistence of the surface 
AP. This functional disconnection between sarcolemma 
and t-system AP generation illustrates that EMG record- 
ings might miss a failure of excitability under conditions 
of high Gm and low /Na(max), as might occur after pro- 
longed muscle activity. 

Regeneration of the t-system AP by voltage-gated ion 
channels was shown in the model to take place at the 
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cost of producing a rapid accumulation of K + ions within 
the t-system during trains of APs. Thus, the K + concen- 
tration within the deepest parts of the t-system [K + ] t ( 19 ) 
was shown to more than double within 2 s of 30-Hz stim- 
ulation. No such K + accumulation was seen in fibers 
lacking voltage-gated Na + channels in the t-system, despite 
the continued presence of voltage-gated K + channels. 
Direct experimental measurement of ion concentra- 
tions within the t-system is difficult or even impossible 
during periods of dynamic change during electrical ac- 
tivity. However, the increased [K + ] t during repetitive AP 
firing was shown to produce a stable depolarization of 
the resting sarcolemmal potential, allowing for compa- 
rison of model findings with experimental microelec- 
trode measurements of the membrane potential. 

Thus, membrane potential measurements were re- 
corded experimentally from rat EDL muscle fibers ex- 
posed to myosin II inhibitor BTS (Macdonald et al., 
2005). The experimental results were in close quanti- 
tative agreement with the model simulations, with 
each showing a baseline resting potential depolariza- 
tion of ~1, 2, and 3 mV during 6-, 15-, and 30-Hz AP 
trains, respectively, and an associated decrease in max- 
imum AP height. The close quantitative agreement 
between model and experiments on the depolariza- 
tion during AP trains permitted the use of the model 
to investigate the underlying alterations in t-system 
K + homeostasis. This revealed that the dynamics of 
t-system K + accumulation during AP trains and recovery 
thereafter are very rapid, completing within 1-2 s. 
This means that measurements of interstitial K + using 
microdialysis probes that offer temporal resolutions in 
minutes might offer rather little information on the 
K + dynamics in the working muscle (Mohr et al., 2004) . 

Combining experiments and simulations further 
demonstrated a significant influence of membrane CP 
permeability upon t-system K + accumulation, membrane 
potential depolarization, and the relationship between 
these two variables. Reduction of extracellular [CI ] by 
replacement with an impermeable anion greatly in- 
creased the magnitude of the depolarization during AP 
trains in both simulations and parallel experiments. 
The GyJ Ga ratio has been measured previously under 
identical experimental conditions (Pedersen et al., 
2009a). Thus, reductions in [CP] e increased depolar- 
ization in both the model and experimental records by 
~50% for 80 mM (Gk/Gq = 0.28) and >100% for 50 mM 
(Gk/ Gci = 0.66) when compared with that in solutions 
of 127 mM {GsJG a = 0.11). Both experimental and 
model muscle fibers showed myotonic behavior in zero- 
[CP] solutions; therefore, a small dose of TTX (10~ 8 M) 
was added in the experiments, and a 50% reduction in 
^Na(max) was applied to the model. Despite this manipu- 
lation, CP-free solutions increased the depolarization by 
410% in the model and 540% in the experimental studies. 
This close agreement between the membrane potential 



records of the model and the experiments prompted 
detailed analysis of the influence of CP reduction upon 
t-system K + homeostasis using the model. Reduced [CP] 
was associated with only minor increases in the concen- 
tration of t-system K + during the AP trains. Thus, under 
normal conditions, CP influences the resting mem- 
brane potential during AP trains by maintaining a high 
Ga/ Gk, such that the influence of t-system K + accumula- 
tion on the surface membrane potential is small. 

Finally, modeling and experimental membrane po- 
tential measurements were used to demonstrate and 
quantify the influence of the physiological Phase I and 
Phase II changes in Ga and Gk upon membrane depo- 
larization and t-system K + accumulation that occur dur- 
ing repetitive AP firing. The Phase I reduction in Ga 
was shown to approximately double the magnitude of 
the resting potential depolarization during 15-Hz stim- 
ulation in both model and experiment, similar to the 
influence of reduced extracellular [CP] . The model 
also indicated an increase in t-system K + accumulation 
from ^10 to ^1 1.5 mM. 

The Phase II increase in Gk and Ga produced a smaller 
increase in baseline depolarization during 15-Hz AP 
trains, although the model slightly underestimated the 
magnitude of the experimentally observed increase. 
The model showed a large reduction in t-system K + accu- 
mulation during Phase II, to a peak of ^7.5 mM com- 
pared with ^10 mM under control conditions. Thus, 
during Phase II, depolarization is increased despite a 
smaller rise in K + . Indeed, a biphasic relationship was 
found for the influence of Ga values on the membrane 
potential during AP firing, such that both decreases 
and increases in Ga in the physiological range increase 
the resting potential depolarization during AP trains. 
The effect of a reduction in Ga was shown to result from 
a shift of the membrane potential closer to which is 
depolarized during activity because of increased [K + ] t . 
In contrast, an increase in G a was shown to produce a 
large contribution of CP currents to AP repolarization, 
thereby causing intracellular CP accumulation and 
hence a depolarization of E a - Phase II also increased 
the rate of recovery to normal values of [K + ] t after activ- 
ity. However, the presence of high background ion con- 
ductances during AP firing in Phase II produced increases 
in both inward and outward K + movements. Thus, al- 
though Phase II Gm changes appear to accelerate reup- 
take of K + from the t-system into the cell and thereby 
enhance recovery of normal ionic concentrations, they 
also appear to reduce the energy efficiency of APs by in- 
creasing the total ionic fluxes required for each AP. 

This study therefore provides a quantitative analysis, 
supported by experimental results, of the relationships 
between resting conductances, excitability, and t-system 
ionic homeostasis in skeletal muscle. It demonstrates 
that active signal regeneration within the t-system pro- 
duces K + accumulation and quantifies the relationship 
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between this and the depolarization of the resting mem- 
brane potential that is observed during trains of APs. It 
shows that the Phase I ^60% decrease in Gm during re- 
petitive AP firing enhances endplate and t-system excit- 
ability, while also increasing the magnitude of t-system 
K + accumulation and approximately doubling the mem- 
brane potential depolarization during AP trains. In con- 
trast, the Phase II increase in Gm decreases endplate 
and t-system excitability and also enhances K + reuptake 
from the t-system, thereby reducing K + accumulation 
during AP trains. Yet, Phase II also increases the resting 
potential depolarization during AP trains because of 
intracellular CI accumulation. 

The model that is described here therefore reproduces 
predictions of electrophysiological behavior and ion con- 
centrations, at rest and during activity, within both the 
sarcoplasm and the t-system, that are in full agreement 
with experimental and prior theoretical results. In so 
doing, it reveals the mechanisms of observed changes in- 
cluding baseline resting potential depolarization during 
AP trains. Furthermore, it provides a platform for further 
exploration of electrophysiological activity in skeletal 
muscle that is applicable where steady-state ion concen- 
trations are not known or expected to undergo signifi- 
cant change. This, in contrast to existing models, permits 
analysis of the influence of activity-related, genetic, and 
pharmacological manipulations to the normal physio- 
logical state. 
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